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ABSTRACT 

We present evidence for a new Milky Way stellar tidal stream in the direction of the Andromeda 
and Triangulum (M31 and M33) galaxies. Using a matched-filter technique, we search the Sloan 
Digital Sky Survey Data Release 8 by creating stellar density maps which probe the Milky Way halo 
at distances between 8 and 40 kpc. A visual search of these maps recovers all of the major known 
stellar streams, as well as a new stream in the direction of M31/M33 which we name the Triangulum 
stream. The stream spans 0.2° by 12° on the sky, or 75 pc by 5.5 kpc in physical units with a best fitting 
distance of 26 ±4 kpc. The width of the stream is consistent with being the tidal remnant of a globular 
cluster. A color magnitude diagram of the stream region shows an overdensity which, if identified as a 
main sequence turn-off, corresponds to an old (~ 12 Gyr) and metal-poor ([Fe/H] ~ — 1.0 dex) stellar 
population. Future kinematic studies of this and similar cold streams will provide tight constraints 
on the shape of the Galactic gravitational potential. 
Subject headings: Galaxy: halo — Galaxy: structure 



1. INTRODUCTION 

Stellar streams are found throughout the galactic 
halos of the Milky Way, Andromeda (hereafter M31) 
and other nearby massive galaxies (e. g., Bclok urqv et al.l 
2001 llbata et al.l [200H: iMartmez-Delgado et al.l I2008D . 
These streams are presumed to originate from disrupting 
dwarf galaxies or globular clusters, providing evidence 
that galaxies are formed, at least in part, hierarchically. 

The recent discovery of numerous Milky Way stel- 
lar streams is due to wide-field photometric surveys 
such as the Sloan Digital Sky Survey (SDSS). Many 
of these streams are associated with known objects 
such as the tidal tail system of the Sagittarius dwarf 
spher oidal galaxy (|Maiewski et a l. 2003, Bclo kurov et al.l 
2006) or with g l obular clusters (e.g., Palomar 5; 
Odenkirchen et all 120031 Palomar 14; S oilima et al.l 



20111) . Some streams, however, are identified as isolated 



features in which the progenitor object is not yet iden- 
tified or has been tidally disrupted beyond re cognition 
(e.g., M onoceros Stream. iNewberg et al~l 120021 : Orphan 
Stream, iBelokurov et al.l 120071 ). 

Stellar streams can extend over large physical dis- 
tances and are sensitive to the strength and shape of 
the Galactic gravita t ional potential (e.g., iHelmil 120041 : 
Uohnston et al.l 120051: ILux et al.ll2012f). as well as to its 
'lumpiness' dCarlberdl2009tlYoon et al.ll201lD . The GD- 
1 tidal stream (jGrillmair fc Dion atos 2006) is sufficiently 
nearby (10 kpc) that proper m o tions are available for its 
member stars. IKoposov et al.l (|2010f ) analyzed the full, 
six-dimensional phase space for GD-1, placing a strong 
constraint on the flattening of the Milky Way potential 
between 10-15 kpc which is dominated by the Galactic 
disk. An analysis of more distant streams will provide 
similar constraints on the Galactic halo potential. Each 
individual stream provides independent constraints on 
the Galactic potential, thus a joint analysis of many 
streams provides a promising method to constrain the 
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full Milky Way potential. 

Here we present evidence for a new Milky Way stellar 
stream, which we name the Triangulum stream, found 
in projection against M31 and M33. In §[2] we discuss 
the SDSS data analysis and matched-filter technique. In 
§[3J we determine the stream properties, and discuss its 
relationship to other Galactic substructures in § 13.31 

2. DATA ANALYSIS 

We search for stellar str eams in the Sloan D igital Sky 
Survey Data Release 8 (|Aihara et al.l 120111 hereafter 
SDSS DR8). We select stars brighter than r = 22 over 
the 14,555 deg 2 SDSS coverage, of which 5200 deg 2 is 
newly released imaging in the Southern Galactic hemi- 
sphere. At these magnitudes, the majority of stars be- 
long to the Milky Way disk. Stars in substructures typi- 
cally belong to different stellar populations as compared 
to the Milky Way disk and can be distinguished by their 
c hemical compo si tion a nd age. 

iRockosi et alj (|2002D first introduced the matched- 
filter technique in the context of Milky Way substruc- 
ture searches. Matched filtering works by preferentially 
weighting stars, in color-magnitude space, with a higher 
probability of being part of a stream population. Fil- 
ters used in this work are based on the globular cluster 
M13, whose old age and low met allicity (13.5 Gyr and 
[Fe/H] = —1.6; iVanden Berg 200(J) are similar to other 
known cold stellar streams. M13 is observed by SDSS, 
thus no filter transformations are required. 

Our filter is constructed by dividing a Hess diagram of 
the target population (Figure [Th.) by a Hess diagram of 
the foreground (FigureO}}). The filter target stars are se- 
lected from the DR8 photometric catalog within a radius 
of 15' from M13's center. The foreground population is 
selected locally, between 45' and 75' from M13. The fi- 
nal filter is smoothed and confined to a region following 
the target population's isochrone within a generous tol- 
erance (Figure [TJ:), which accommodates isochrones of 
other commonly used globular clusters (e.g., M92). 

Our approach ensures that regions in the color- 
magnitude diagram (CMD) with high contrast between 



the target and foreground population are given the most 
weight in filtering. The filter is most heavily weighted 
near the main sequence turn-off (Figure [TJ:) , where nu- 
merous Mf 3 stars are present, but few Galactic stars are 
as blue in color. The main sequence is less weighted be- 
cause of the large population of Milky Way stars present 
redward of the turn-off. The subgiant and giant branches 
are weighted the least because of the small number of 
M13 stars which occupy those areas of the CMD. 

A total of four matched filters (in r vs. g — r, g vs. g — r, 
g vs. g — i and g vs. g — z) were applied to the whole 
SDSS survey area for a distance range between 8 kpc and 
40kpc. The inner distance limit was chosen to avoid the 
Milky Way disk, while the outer limit is the distance at 
which the main sequenc e turnoff falls bel ow the SDSS 
95% completeness limit (|Juric et all 12008( 1 . Each filter 
was applied to the full SDSS photometric sample and the 
resulting density maps were coadded. The left column 
of Figure [2] shows the final density maps for the North 
Galactic hemisphere, while right panels show the South 
Galactic hemisphercQ The SDSS footprint in the North 
has been extensively searched for halo substructure and i s 
known as the "field of streams" (IBelokurov et al.ll2006f l. 
We recover both compact and diffuse known objects in 
this field. Globular clusters and dwarf galaxies are visible 
as small, concentrated overdensities. At a distance of 
8 kpc (Figure [2J top) we detect the cold GD-f stream, 
as well as the much wider signatures of the Monoceros 
stream and leading Sagittarius tidal arm. At a distance 
of 26 kpc (Figure EJ bottom) we detect the tidal tails of 
the Palomar 5 globular cluster and the Orphan stream. 
A visual search of these maps does not reveal any new 
stellar streams beyond those which have been published. 
Our method recovers streams down to the li mit of the 
faint Lethe and Styx streams (|Grillmairll2009[ ). 

Stellar density maps of the South Galactic hemisphere 
are shown in the right column of Figure [5J at distances of 
8 kpc (top) and 26 kpc (bottom) . The Sagittarius trailing 
tidal arm is visible at both distances. At larger latitudes 
we also detect a fainter strea m parallel to th e stron ger 
Sagittarius tail, as discussed in lKoposov et al.l ()2012| ). A 
visual search of the Southern maps reveals a new stream 
in the region of M3f/M33 which we describe in detail 
below. No other substructure features were detected. 

3. RESULTS 

We demonstrate in Figure [2] that our search method 
recovers all known stellar streams thus far published. In 
addition, we find evidence for a new narrow stream in 
the direction of the M31/M33 region which we name the 
Triangulum stream. In § 13.31 we discuss this new stream 
in the context of other known substructure in this region. 

3.1. Physical Properties of the Triangulum Stream 

The Triangulum stream is long and very narrow (Fig- 
ure [3^,). To ensure that the feature is indeed a stel- 
lar overdensity, we s how the dust extinction (Figure [3Jd; 
iSchlegel et al.l fl998| l and number counts of blue SDSS 
galaxies (Figure Et) in the stream region. Neither dust 
nor galaxies show an overdensity in the stream region. 

2 Our S DSS DR8 stellar density m aps are ava ilable at 
http: / /www. astro, yale.edu/abonaca/research/halo. html 



TABLE 1 
Triangulum stream properties 



Property Value 

RA (deg^_ [21.346, 23.975] 

Dec (deg^ [34.977, 23.201] 

I (deg) [130.756, 136.023] 

b (deg) [-27.378, -38.538] 

Width 0.2° (75 pc) 

Length 12° (5.5 kpc) 

Distance 26 zt 4 kpc 

[Fe/H] -1.0 dcx 

Age 12 Gyr 



a Stream boundary points, see Equation [T] for a relation between 
RA and Dec. 

The Triangulum stream is nearly linear on the sky, 
whose coordinates satisfy the equation: 

Dec(°) = -4.4 x RA(°) + 128.5 (1) 

in the range RA~ 21° - 24° (Table [TJ). We define a new 
coordinate system such that the stream runs along a y- 
axis. This system is an equatorial system rotated 11° 
counterclockwise around RA = 22° and Dec = 32°. A 
stellar density map of the stream region in this new co- 
ordinate system is shown in Figure [3ji. To assess the 
statistical significance of this overdensity, we create a 
one-dimensional stellar density profile of the stream (re- 
gions —9° < y < —6° and 0° < y < 3°, marked with 
red and blue lines on Figure [5Ji, respectively) extending 
\x\ < 5° on each side of the stream. For x < —1° the 
stream profile is affected by the incomplete SDSS cover- 
age. We correct for completeness by dividing the profile 
value with the fraction of populated pixels in each ver- 
tical bin. This procedure enables easier analysis of the 
stream profile, but it also introduces additional noise to 
the profile at x < —1°. The resulting profile is shown 
in Figure [3^: blue and red lines are profiles of top and 
bottom stream regions, respectively, while the black line 
is their average. 

The stream signal was strongest for the matched filter 
at 26 kpc. We show this quantitatively by fitting a Gaus- 
sian function with a constant term to the stream profile 
at a range of distances between 20 kpc and 30 kpc. A 
Levenberg-Marquardt minimization algorithm was used 
to solve for the mean, dispersion, and normalization of 
the distribution. Both the dispersion and peak of the 
profile are maximized at a distance of ~ 26 kpc, while the 
mean does not depend on the distance. The height of the 
stream peak is 3 a above the background noise (estimated 
in the region 0.5° < \x\ < 5°) at distances 23 — 30 kpc. 
We adopt a distance to the stream of 26 ± 4 kpc. The 
measured angular extent is A x = 10' and A y = 12°, 
which correspond to a width of 75 pc and a length of 
5.5 kpc in physical units. 

We also detected a ~ 2 a overdensity of blue horizon- 
tal branch (BHB) stars in the lower part of the stream 
region (-9° < y < -6°). In this box, the 0.5° wide 
strip contains on average 7.5 ± 0.7 BHBs in the stream 
fields (\x\ < 0.5°) and 4.6 ± 2.1 BHBs in the off fields 
(0.5° < \x\ < 5°). 

In general, the shape of a stellar stream traces the 
past orbit of its progenitor. The stream width corre- 
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Fig. 1. — Steps in creating an M13-based filter for selection of old, metal-poor stars: the left panel shows a Hess diagram of SDSS 
stars detected in the inner region of globular cluster M13, while that of a surrounding area, representative of Milky Way foreground, is 
given in the middle panel. A ratio of the two produces the desired filter, which preferentially selects stars that are in the region of the 
color-magnitude diagram with the greatest contrast between M13 and Milky Way stellar populations. In order to reduce dumpiness, the 
filter is smoothed (right panel). The resulting filter is most heavily weighted near the main sequence turn-off. Red lines, derived from an 
isochrone matching the M13 population, trace the width of the filter on all three panels. 
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Fig. 2. — Stellar density maps obtained via matched filtering at distances of 8kpc (top) and 26kpc (bottom). The left columns show 
the SDSS footprint in the northern Galactic hemisphere, while maps on the right are in the south. The stretch is linear, with darker areas 
corresponding to higher stellar density regions. The most prominent features are tidal arms of the Sagittarius dwarf galaxy, bright Milky 
Way globular clusters, and brightest neighboring galaxies. Regions not imaged by the SDSS are shown in white. The maps are available 
online. 



lates with progenitor mass, while the length of a stream 
is related to the time sinc e the stars became unbound 
(| Johnston fc Bullock] l200l . The average width of the 
Triangulum stream, a — 10' (75 pc), is comparable to 
other globular cluster streams in the literature. The GD- 
1 stream, whose progenitor is undetecte d but presumed 
to be a globular cluster, is a — 12' (35 pc: IKoposqv et al.1 
20101) . while the globular cluster Pal 5 is a = 18' (120 pc; 



Odenkirchen et all 120031 ). In contrast, dwarf galaxies 



produce much broade r streams : the Orphan stream 
is ~ 2° wide (730 pc; iGrillmairl [20061 IBelokurov et~aTl 
2006), while the Sagittarius stream is 20° or more (7kpc; 



IBelokurov eFalll^fM iMaiewski et al.ll2003D . The width 
of the Triangulum stream suggests that its progenitor is 
a globular cluster, consistent with the age and metallicity 
determined below. 

Our estimated length for the Triangulum stream, 12° 
(5.5 kpc), is likely an underestimate given the incom- 
pleteness of the SDSS data. This is shorter, but compara- 
ble to the Palomar 5 stream (22°, 9 kpc). Deeper imaging 
along the Triangulum stream may reveal a longer extent 
than detected in the SDSS data. 

3.2. Age and Metallicity of the Triangulum Stream 



3 




RA (deg) RA (deg) RA (deg) x (deg) 



Fig. 3. — (a) Matched filter map of the Triangulum stream region in equatorial coordinates, (b) Dust map of the same region. Color 
scale represents reddening E(B — V). (c) Blue galaxy number counts in the stream region, (d) Matched filter map of the wider stream 
region in rotated coordinates. Stream regions are boxed in blue and red, the dashed ellipse shows the extent of RGB stars around the M33 
disk. Red lines on panels a-d bracket stream region, (e) ID profiles of the stream in the top and bottom boxes from the panel (d) (blue 
and red, respectively) and their average (black). 



We estimate the age and metallicity of the Triangulum 
stream using the best-fitting distance of 26 ± 4kpc. Due 
to the small number of stars associated with the stream, 
the CMD in this region is dominated by Milky Way 
foregr ound stars. We therefore follow IKoposov et al.l 
(2010) in constructing a statistically subtracted CMD. 
First, we create a ID profile of all stars satisfying the 
color- magnitude cut 20 < r < 22, 0.2 < g — r < 0.5 
in the stream region (|x| < 5°, —9° < y < —6° and 
0° < y < 3°). This produces a well populated profile, 
with predominantly stream stars. Next, we fit a Gaussian 
to this profile to determine the central position xq and 
width (To- To construct a statistically-subtracted Hess 
diagram, we create a ID profile for stars in every pixel in 
the Hess diagram, fit the profile with a Gaussian of cen- 
tral position xg and width oo, and assign each pixel the 
fitted Gaussian normalization. Thus, higher-value pix- 
els correspond to regions with higher density profile in 
the stream region and, consequently, higher probability 
of being occupied by a stream rather than a foreground 
star. 

The statistically subtracted r vs. g — r Hess diagram is 
noisy, but clearly shows an overdensity at r > 20.5 and 
g — r< 0.5 (Figure[4ji). At a distance of 26 kpc, this over- 
density corresponds to main sequence stars, and a hint 
of a main sequence turnoff is indeed seen. The contour of 
a smoothed Hess diagr am is overplotted w ith a range of 
Dartmouth isochrones (jDotter et al.l l2~008. red isochrones 
on Figure Hb). The isochrones are 12 Gyr old and at a 
distance of 26 kpc. The left-most isochrone is the most 
metal-poor, [Fe/H] = —1.5, while the right-most is the 
most metal rich, [Fe/H] = —0.5. These values bracket 
the likely stream metallicity range, with the most proba- 
ble value being the median [Fe/H] = —1.0. If a range of 
13 Gyr isochrones is used, the inferred metallicity is lower 



[Fe/H] = —1.25. A range of younger and more metal rich 
isochrones fails to match the shape of the overdensity in 
the Hess diagram. 

Given the proximity of the stream to the M33 
galaxy, we have also overplotted a range of old and 
metal -poor isochrones a t the distance of M33 (809 
kpc; rMcConnachi e et"ail [2005) to test whether the 
stream could be associated with the M33 system (gray 
isochrones on Figure Hfc). The tip of the red giant branch 
of that isochrone only reaches r = 21.5 at g — r ~ 1, 
much redder than the stream feature. We conclude that 
the stream is not associated with M33, but is most likely 
the tidal remnant of an old, metal-poor globular cluster 
in the Milky Way halo. 

3.3. Relationship to Other Substructures 

The Triangulum stream is in a region of sky rich with 
substructure. While the stream is projected against the 
M31/M33 system, it has no spatial overlap with any 
of the M31 satellites ([Richardson et al.ll2011| ) or diffuse 
structures discovered around M33 (marked with a dashed 
ellipse in FigureEJi). The M33 population is mo re metal- 
poor ([Fe/H] = -1.6; lMcConnachie et alll2010| ) than the 
Triangulum stream. Thus, we believe that Triangulum 
stream is a Galactic feature, unassociated with the sub- 
structure in the M31/M33 system. 

The largest Milky Way feature in this region is the dif- 
fuse Triangulum-Andromeda (TriAnd) overdensity cov- 
ering - 1000 deg 2 between 10 0° < I < 150°, -20 ° > b > 
-40°. It was discovered by IRocha-Pinto et al.l (|2004T ) 
as a smooth overdensity of M-giants in the 2MA SS sur- 
vey between 15 — 30 kpc. iMaiewski et al.l (|2004D refined 
this distance estimate using main sequence turnoff stars 
to 16 — 25 kpc. The mean spectroscopic metallicity of 
TriAnd stars is [Fe/H] = -0.64 ± 0.08, with a CMD in- 
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Fig. 4. — Top: Hess diagram of the stream region with statisti- 
cally subtracted background. Bottom: Contours of the smoothed 
Hess diagram with overplotted 12 Gyr isochrones spanning a metal- 
licity range of [Fe/H] = —0.5 to —1.5 at the stream distance (red) 
and [Fc/H] = -1.1 to -2.1 at the M33 distance (gray). 



ferred age of 8 Gyr (|Chou et all 120111) . The metallicity 
of TriAnd is higher than that inferred for the Triangu- 
lum st ream and we estima te a population which is 4 Gyr 
older. iMartin et al.l (|2007D have detected a signature of 
another overdensity behind the TriAnd, at heliocentric 
distance ~ 28 kpc. Both features are far more extended 
than Triangulum, however, the Triangulum stream could 
be related to these ovcrdensities as a localized feature at 
the outer edge of a much larger diffuse merger remnant. 
Spectroscopic and/or astrometric follow-up is required to 
explore this possibility. 

Recently, it has been argued that the 'bifu rcation' in 
the Sagittarius stream (Bclokur ov et al.l 120061 ) is in fact 
a di stinct stream with it s own progenitor and kinemat- 
ics (|Koposov et al.l I2012D . The so-calle d Cetus stream, 
first discovered bv lNewberg et al.l (|2009f ). is fainter, thin- 
ner and more metal-poor than the Sagittarius stream. 
Its derived orbit is in the direction of the Triangulum 
stream, however, it has lower metallicity ([Fe/H] = —2.1) 
and appears to be at larger heliocentric distance (34 kpc; 
iNewberg etldl I2009D . Similarly, we do not think that 
the Triangulum Stream is related to the Monoceros 
Ring/ Anticenter Stellar Structure/Monoceros Stream, 
which is the other prominent diffuse structure in the 



South. Even though the latter structure (whose origin 
remains controversial) has a large metallicity spread from 
[Fe/H] = -1.6 to [Fe/H] = -0.4 that could be consistent 
with the stream discovered here, it seems t o be at much 
close r heliocentric distance at ~ 11 kpc (jYannv et al.l 
llMlcTane et al.ll200l IChou et al.ll2010l) . 

4. CONCLUSIONS 

We have searched the SSDS DR8 photometric database 
for Milky Way stellar streams using a matched filtering 
technique. We recover all known streams in the current 
literature. Based on a visual search of our stellar density 
maps, we present evidence for a new stream discovered 
in the southern Galactic hemisphere in the direction of 
M33 which we name the Triangulum stream. We did not 
find any other stream-like structures in the dataset. 

The Triangulum stream spans 0.2° by 12° (75 pc by 
5.5 kpc) on the sky with a best fitting distance of 
26 ± 4 kpc. We estimate the age and metallicity of the 
stream via a statistically subtracted color-magnitude di- 
agram of the region, which suggests a 12 Gyr, metal-poor 
([Fe/H] = —1.0) stellar population. The narrow width 
of the Triangulum stream is comparable to the GD-1 
and Palomar 5 streams, from which we conclude that 
the Triangulum stream is also a disrupted globular clus- 
ter. While we do not observe the progenitor object, the 
stream extends beyond the SDSS ph otometric coverage. 
We n ote that the PAndAS team dMcConnachi e et al.l 
2009) has independently detected the Triangulum stream 
at the same sky position but at a higher significance due 
to deeper photometric data (M. Fardal, priv. communi- 
cation) . 

Searches for cold stellar streams in the halo are moti- 
vated by the desire to map the Galactic gravitational po- 
tential. Further characterizing the size and shape of the 
Triangulum stream via deeper and more extended imag- 
ing will improve the usefulness of this stream as a Galac- 
tic halo probe. Radial velocities and proper motions of 
Triangulum stream stars will provide tighter constraints 
on the nature of its stellar populations and orbital prop- 
erties, providing new constraints on the Galactic gravi- 
tational potential. 
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